Best practices for genotype imputation from low-coverage sequencing data in natural populations.

Monitoring genetic diversity in wild populations is a central goal of ecological and evolutionary genetics and is critical for conservation biology. However, genetic studies of nonmodel organisms generally lack access to species-specific genotyping methods (e.g. array-based genotyping) and must instead use sequencing-based approaches. Although costs are decreasing, high-coverage whole-genome sequencing (WGS), which produces the highest confidence genotypes, remains expensive. More economical reduced representation sequencing approaches fail to capture much of the genome, which can hinder downstream inference. Low-coverage WGS combined with imputation using a high-confidence reference panel is a cost-effective alternative, but the accuracy of genotyping using low-coverage WGS and imputation in nonmodel populations is still largely uncharacterized. Here, we empirically tested the accuracy of low-coverage sequencing (0.1-10×) and imputation in two natural populations, one with a large (n = 741) reference panel, rhesus macaques (Macaca mulatta), and one with a smaller (n = 68) reference panel, gelada monkeys (Theropithecus gelada). Using samples sequenced to coverage as low as 0.5×, we could impute genotypes at >95% of the sites in the reference panel with high accuracy (median r2 ≥ 0.92). We show that low-coverage imputed genotypes can reliably calculate genetic relatedness and population structure. Based on these data, we also provide best practices and recommendations for researchers who wish to deploy this approach in other populations, with all code available on GitHub (https://github.com/mwatowich/LoCSI-for-non-model-species). Our results endorse accurate and effective genotype imputation from low-coverage sequencing, enabling the cost-effective generation of population-scale genetic datasets necessary for tackling many pressing challenges of wildlife conservation.

Vascularization underlies differences in sexually selected skin coloration in a wild primate.

Male reproductive competition can select for condition-dependent, conspicuous traits that signal some aspect of fighting ability and facilitate assessment of potential rivals. However, the underlying mechanisms that link the signal to a male's current condition are difficult to investigate in wild populations, often requiring invasive experimental manipulation. Here, we use digital photographs and chest skin samples to investigate the mechanisms of a visual signal used in male competition in a wild primate, the red chest patch in geladas (Theropithecus gelada). We analysed photographs collected during natural (n = 144) and anaesthetized conditions (n = 38) to understand variability in male and female chest redness, and we used chest skin biopsies (n = 38) to explore sex differences in gene expression. Male and female geladas showed similar average redness, but males exhibited a wider within-individual range in redness under natural conditions. These sex differences were also reflected at the molecular level, with 10.5% of genes exhibiting significant sex differences in expression. Subadult males exhibited intermediate gene expression patterns between adult males and females, pointing to mechanisms underlying the development of the red chest patch. We found that genes more highly expressed in males were associated with blood vessel development and maintenance but not with androgen or oestrogen activity. Together, our results suggest male gelada redness variability is driven by increased blood vessel branching in the chest skin, providing a potential link between male chest redness and current condition as increased blood circulation to exposed skin could lead to heat loss in the cold, high-altitude environment of geladas.

Urinary neopterin reflects immunological variation associated with age, helminth parasitism, and the microbiome in a wild primate.

Neopterin, a product of activated white blood cells, is a marker of nonspecific inflammation that can capture variation in immune investment or disease-related immune activity and can be collected noninvasively in urine. Mounting studies in wildlife point to lifetime patterns in neopterin related to immune development, aging, and certain diseases, but rarely are studies able to assess whether neopterin can capture multiple concurrent dimensions of health and disease in a single system. We assessed the relationship between urinary neopterin stored on filter paper and multiple metrics of health and disease in wild geladas (Theropithecus gelada), primates endemic to the Ethiopian highlands. We tested whether neopterin captures age-related variation in inflammation arising from developing immunity in infancy and chronic inflammation in old age, inflammation related to intramuscular tapeworm infection, helminth-induced anti-inflammatory immunomodulation, and perturbations in the gastrointestinal microbiome. We found that neopterin had a U-shaped relationship with age, no association with larval tapeworm infection, a negative relationship with metrics related to gastrointestinal helminth infection, and a negative relationship with microbial diversity. Together with growing research on neopterin and specific diseases, our results demonstrate that urinary neopterin can be a powerful tool for assessing multiple dimensions of health and disease in wildlife.

A novel statovirus identified in fecal samples from wild geladas in the Ethiopian highlands.

We present a novel statovirus in geladas (Theropithecus gelada), graminivorous primates endemic to the Ethiopian highlands. Using a high-throughput sequencing approach, we identified contiguous sequences in feces from two adult female geladas in the Simien Mountains National Park, Ethiopia, that share similarities to statoviruses. Our phylogenetic analysis of the whole genome, as well as the RNA-dependent RNA polymerase (RdRp) and capsid protein (CP) amino acid sequences, revealed that the gelada statoviruses cluster with those from other primates (laboratory populations of Macaca nemestrina and Macaca mulatta). As the first report of statovirus in wild primates, this finding contributes to our understanding of the phylogenetic and geographic distribution of statoviruses and their hosts.

Genomic signatures of high-altitude adaptation and chromosomal polymorphism in geladas.

Primates have adapted to numerous environments and lifestyles but very few species are native to high elevations. Here we investigated high-altitude adaptations in the gelada (Theropithecus gelada), a monkey endemic to the Ethiopian Plateau. We examined genome-wide variation in conjunction with measurements of haematological and morphological traits. Our new gelada reference genome is highly intact and assembled at chromosome-length levels. Unexpectedly, we identified a chromosomal polymorphism in geladas that could potentially contribute to reproductive barriers between populations. Compared with baboons at low altitude, we found that high-altitude geladas exhibit significantly expanded chest circumferences, potentially allowing for greater lung surface area for increased oxygen diffusion. We identified gelada-specific amino acid substitutions in the alpha-chain subunit of adult haemoglobin but found that gelada haemoglobin does not exhibit markedly altered oxygenation properties compared with lowland primates. We also found that geladas at high altitude do not exhibit elevated blood haemoglobin concentrations, in contrast to the normal acclimatization response to hypoxia in lowland primates. The absence of altitude-related polycythaemia suggests that geladas are able to sustain adequate tissue-oxygen delivery despite environmental hypoxia. Finally, we identified numerous genes and genomic regions exhibiting accelerated rates of evolution, as well as gene families exhibiting expansions in the gelada lineage, potentially reflecting altitude-related selection. Our findings lend insight into putative mechanisms of high-altitude adaptation while suggesting promising avenues for functional hypoxia research.

Helminth infection is associated with dampened cytokine responses to viral and bacterial stimulations in Tsimane forager-horticulturalists.

BACKGROUND: Soil-transmitted helminths (STHs) and humans share long co-evolutionary histories over which STHs have evolved strategies to permit their persistence by downregulating host immunity. Understanding the interactions between STHs and other pathogens can inform our understanding of human evolution and contemporary disease patterns. METHODOLOGY: We worked with Tsimane forager-horticulturalists in the Bolivian Amazon, where STHs are prevalent. We tested whether STHs and eosinophil levels-likely indicative of infection in this population-are associated with dampened immune responses to in vitro stimulation with H1N1 and lipopolysaccharide (LPS) antigens. Whole blood samples (n = 179) were treated with H1N1 vaccine and LPS and assayed for 13 cytokines (INF-γ, IL-1ß, IL-2, IL-4, IL-5, IL-6, IL-7, IL-8, IL-10, IL-12p70, IL-13, GM-CSF and TNF-ɑ). We evaluated how STHs and eosinophil levels affected cytokine responses and T helper (Th) 1 and Th2-cytokine suite responses to stimulation. RESULTS: Infection with Ascaris lumbricoides was significantly (P ≤ 0.05) associated with lower response of some cytokines to H1N1 and LPS in women. Eosinophils were significantly negatively associated with some cytokine responses to H1N1 and LPS, with the strongest effects in women, and associated with a reduced Th1- and Th2-cytokine response to H1N1 and LPS in women and men. CONCLUSIONS AND IMPLICATIONS: Consistent with the 'old friends' and hygiene hypotheses, we find that STHs were associated with dampened cytokine responses to certain viral and bacterial antigens. This suggests that STH infections may play an essential role in immune response regulation and that the lack of STH immune priming in industrialized populations may increase the risk of over-reactive immunity. Lay Summary: Indicators of helminth infection were associated with dampened cytokine immune responses to in vitro stimulation with viral and bacterial antigens in Tsimane forager-horticulturalists in the Bolivian Amazon, consistent with the 'old friends' and hygiene hypotheses.

Molecular identification of Taenia serialis coenurosis in a wild Ethiopian gelada (Theropithecus gelada).

Since morphological identification of a larval Taeniid in geladas (Theropithecus gelada) has produced inconsistent results, genetic information is pivotal for species identification. Nuclear and mitochondrial DNA from a coenurus in a wild gelada were compared to published sequences from multiple Taeniid species, confirming the identification of this parasite as Taenia serialis. A demographic analysis finds age to be a strong predictor of coenuri. Tapeworms rarely employ primates as intermediate hosts, and the presence of T. serialis in a wild gelada population may indicate a substantial ecological shift in this parasite's life cycle.